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ABSTRACT 

Reverberation lags in AGN were first discovered in the NLSl galaxy, 1H0707-495. We 
present a follow-up analysis using 1.3 Ms of data, which allows for the closest ever look 
at the reverberation signature of this remarkable source. We confirm previous findings 
of a hard lag of 100 seconds at frequencies v ^ [0.5 — 4] x 10^"* Hz, and a soft lag 
of ^ 30 seconds at higher frequencies, v ~ [0.6 — 3] x 10^"^ Hz. These two frequency 
domains clearly show different energy dependences in their lag spectra. We also find 
evidence for a signature from the broad Fe Ka line in the high frequency lag spectrum. 
We use Monte Carlo simulations to show how the lag and coherence measurements 
respond to the addition of Poisson noise and to dilution by other components. With our 
better understanding of these effects on the lag, we show that the lag-energy spectra 
can be modelled with a scenario in which low frequency hard lags are produced by 
a compact corona responding to accretion rate fluctuations propagating through an 
optically thick accretion disc, and the high frequency soft lags are produced by short 
light-travel delay associated with reflection of coronal power-law photons off the disc. 

Key words: black hole physics - galaxies: active - X-rays: galaxies - galaxy: indi- 
vidual : 1H0707-495. 



X: 
S3: 



1 INTRODUCTION 

Recent X-ray variability studies of active galactic nuclei 
(AGN) have given us unique insights into the energetics 
and geometry of the innermost regions of accreting black 
hole systems. By measuring the phase lags between light 
curves at different energy bands, we can probe distances as 
small as tens of lightseconds in AGN that are billions of 
lightyears away. Narrow-line Seyfert I galaxy 1H0707-495 
{z = 0.0411) was the first AGN shown to exhibit rever- 
beration lags (jPabian et al.l |2009| ). and since then, this 
feature has bee n significantly detected in 15 other Seyfert 
galax ies temm anoulopoulos. M cHardv. fc PapadakisI 



20 111: Ide Mar co ct alf 



2011 



Zoghbi & F abianI l2oTir 



de Marco et all l2012i : IZoghbi et aLri2012l ). Further 



dence for reverberation lags of milliseconds or less have 
been seen in stellar mass black hole X-ray binary GX 339-4 
(jUttlev et al.ll201ll '). 

Reverberation lags can be explained by the standard 
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reflection model. Matter in a geometrically thin, optically 
thick accretion disc is accreted onto the black hole, caus- 
ing the accreti on disc to heat up and ra diate, mostly in 
the ultraviolet (|Shakura fc SunvaevI Il973 ). Those thermal 
disc photons are Compton upscattered to X-ray energies 
in the corona or base of the jet, producing the observed 
X-ray power-law continuum. Some continuum photons are 
directed back towards the accretion disc, causing it to fluo- 
resce with the characteristic reflection spectrum. Soft lags or 
'reverberation lags', where the power-law continuum leads 
the 'soft excess' associated with the photoionized disc re- 
flector, are thought to be due to the effect of the light 
travel time between the corona and the inner accretion disc. 
Hard lags, where the power-law continuum follows varia- 
tions at soft energies, have been explained as the propaga- 
tion of mass accretion rate fluctuations, which modulate re- 
gions emitting soft photons before those emitting hard pho- 
tons jKotov. Churazov. fc Gilfanov|[200ll : lArevalo fc UttlevI 
l2006h . 

XMM-Newton observations of 1H0707-495 have been 
some of the most important resources in gaining an under- 
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standing of the X-ray emission from the innermost regions 
of active galaxies. The first XMM-Newton observations of 
1H0707-495 showed a sharp drop in the spectrum at ~ 7 keV 
l|Boller et al.ll2002l V Two interpretations were developed to 
explain this feature, n amely an absorption edge from a par- 
tial covering absorber (iTanaka et al.ll2004l ') , and the edge to 
the blue wing of a broad iron Kg line, produced in the inner 
regions of the disc (jPabian et al.|[20M ). 

Further observations in 2008 showed more evidence for 
the inner disc interpretation, as the co mplementary broad 
iron La line was a lso clearly identified ijFabian et al.|[2009l : 
IZoghbi et alj [2O10h . In addition to this spectral evidence, 
timing studies showed a ~ 30 second light travel delay 
betwee n the pow er law source and the reflection from the 
disc. .Miller et al.l (2010) interpreted this time delay as be- 
ing due to scattering in a more distant reflector/absorber 
associated with a wind directed along the line of sight, 
though the subsequent X-ray reverberation detections in 
several different AGN make this an unlikely scenario. See 
IZoghbi. Uttlev. fc Fabia n (2011) for further discussion re- 
garding problems with this model. 

In 2011, 1110707-495 went into a low flux state, which 
triggered an observation with XMM-Newton. The spectral 
drop associated with the blue wing of the Fe Ka appeared 
at ~ 6 keV, suggesting emission very close t o the central 
black hole, within 1 rg of the event horizon (iFabian et al.l 
|2012| ). Lastly, recently public observations from 2010 confirm 
refiection from a highly ionized disc, and suggest that the 
reflection spectr um could b e slightly modified by a highly 
ionized outfiow (jPauser et al. 2012). 

In this paper, we compile all the XMM-Newton observa- 
tions from 2000 to 2010 (> 1.3 Ms), to take the closest ever 
look at the reverberation lags in 1H0707-495. With ~ 2.5 
times the dataset, we are better able to define the lag spec- 
trum and s tructures in the data that were s uggested by the 
analysis of IZoghbi. Uttlev. fc FabianI l|201ll ). The observa- 
tions and data reduction are described in Section [51 and 
the results are presented in Section [3] We continue with a 
study of the robustness of the lag measurement using Monte 
Carlo simulations in Section 3] and discuss our interpreta- 
tion through a physical model for the lag spectrum in Section 
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Table 1. The 14 JfMM-Afeuiton observations used in this analysis, 
amounting to 1.3 Ms of data. 

node as the source region. In order to get the best descrip- 
tion of the background, we chose the background regions to 
be as large as possible, often twice the radius of the source 
region. The position of the background regions were chosen 
to avoid the Cu-K emission lines from the electronic circuits 
behind the PN CCD that contaminate the background at 8.0 
and 8.9 keV. Despite these efforts, the Cu-K line still mildly 
contaminates the data, and so we take heed with regards to 
the highest energy bin at 7-10 keV. 

The response matrices were produced using rmfgen 
and ARFGEN in SAS. The spectra from all observations were 
merged before fitting using mathpha in FTOOLS, and the 
resulting combined spectra were rebinned to contain a min- 
imum of 20 counts per bin. Because the observations differ 
significantly in flux, we also used tabgtigen in SAS to filter 
the observations by rate, in order to check the flux resolved 
spectra and compare them to the combined spectrum from 
the 14 observa tions. Spectra l fitting was performed using 
XSPEC vl2.5.0 (|Arnaudll 19961 ). All quoted errors correspond 
to a 90% confidence level, and energies are given in the rest 
frame of the so urce. The quo ted abundances refer to solar 
abundances in I Anders fc GrevessQ (,19891 '). 



2 OBSERVATIONS AND DATA REDUCTION 

1H0707-495 has been observed several ti mes over the past 1 2 
years with the XMM-Newton satellite (|jansen et al.ll200lf ). 
The 14 observations used in this analysis are shown in Ta- 
ble [J We focus on the data from the EPIC-PN camera 
l|Striider et al.ll200ll ). which were all made in large window 
imaging mode, with the exception of the 2000 and 2001 data, 
which were taken in full window imaging mode. All the data 
were reduced in the same way using the XMM-Newton Sci- 
ence Analysis System (SAS v. 11. 0.0) and the newest cal- 
ibration files. The data were cleaned for high background 
fiares, and were selected using the condition pattern ^ 4. 
No significant pile-up was found in any of the observations. 

The source spectra were extracted from circular regions 
of radius 35 arcsec centered on the maximum source emis- 
sion, and the background spectra were chosen from back- 
ground regions that were the same distance to the readout 



3 RESULTS 

3.1 Lag vs. Frequency Spectrum 

Using the light curves of the 14 observations, we compute 
the Fourier phase lag between the hard an d soft energy 
bands, following the technique described in iNowak et al.l 
(1999). We write the Fourier transform of the soft light 
curve, s, as 5* = |5'|e*'''''. The product of S with the com- 
plex conjugate of the hard band FFT, H* , is written, 
H'S = |H||S'|e'('*=-'*''', which gives us the phase difference 
between the two signals. The overall Fourier phase lag, (j){f) 
is the phase of the average cross power spectrum. That is, 
(f>{f) — arg[(//*(/)5'(/))]. We then convert this to a time 
lag, T = (j>{f) /"ii^f- We note that the Fourier phase is defined 
on the interval ( — tt, tt), which can cause p hase wrapping at 
high frequencies l|Nowak fc Vaughanlll99(j '). 

Fig. [T] shows the lag-frequency spectrum between the 
energy bands 0.3 - 1.0 keV (soft-excess dominated) and 



© 2012 RAS, MNRAS 000, [HIT] 



The Closest Look at 1H0707-495: X-ray Reverberation Lags with 1.3 Ms of Data 3 



200 




0.01 



Fourier Frequency (Hz) 



Figure 1. Lag-frequency spectrum using 1.3 Ms of data. The lag 
is calculated between the soft energy band (0.3 - 1. keV) and the 
hard band (1.2 - 4. keV). We adopt the convention that negative 
lag mean the soft band lags behind the hard band. The most 
negative lag (at 1.33 X 10"^ Hz) is -31.9 ± 4.2 s. 



1.2 - 4.0 keV (power law dominated). The light curve seg- 
ments range in individual exposure times from 40000 sec- 
onds to 110000 seconds with 10 second bins. The frequency 
range goes up to 0.01 Hz, beyond which the signal is dom- 
inated by Poisson noise. As seen in p re vious timing s tudies 
of 1 H0707-495 bviFabian et all ll2009h. IZoehbi et all (l2010l ) 



and lEmmanoulopoulos. McHardv. fc Papadakid 1 2011 ). the 
hard flux is shown to lag behind the soft by hundreds of 
seconds at frequencies less than ~ 5 x lO"** Hz (here- 
after referred to as the "hard lag"). Above frequencies of 
~ 7 X lO"** Hz, the soft flux is clearly shown to lag the hard 
by roughly 30 seconds (referred to as the "soft lag"). 

3.2 Lag vs. Energy Spectrum 

Taking a closer look, we examine how the time lag evolves 
with energy, using roughly equal logarithmically-spaced en- 
ergy bands (henceforth 'ener gy bins') in a particul ar fre - 
quency range, similar to iZoghbi. Uttlev. fc Fabiani (|201ll ). 
For each energy bin, the lag is computed between the light 
curve in that bin and the light curve of a broader refer- 
ence band from 0.3 - 0.8 keV. This reference energy range 
is selected because it is the band with the highest signal- 
to-noise. For lags computed in energy bins between 0.3 - 
0.8 keV, the light curve of the energy bin is removed from 
the reference to ensure that the Poisson noise remains un- 
correlated. The lags have not been shifted, and therefore 
zero-lag means that there is no time delay between that bin 
and the soft excess reference band. Similarly, a negative lag 
means that that bin leads the reference, and positive lags 
mean that that bin la gs behind the referenc e . Prev ious lag- 
energy studies by .Zoghbi. Uttlev. fc Fabiani (|201ll ) use the 
entire energy range as the reference band. We find that this 
does not change the shape of the lag-energy profile, but does 
make for an arbitrary constant lag shift in the entire spec- 
trum. By choosing the reference band from 0.3 - 0.8 keV, 
we have a more physical interpretation of what the value of 
the lag means, namely that it lags or leads the soft excess. 
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Figure 2. Lag-energy spectra showing the energy dependence of 
the lags at particular frequencies, (a) The lag-energy spectrum for 
very-low frequencies, u = [0.5 — 1.65] X lO"** Hz (blue circle; left 
axis), and low frequencies, u = [1.65 — 5.44] X 10~* Hz (red square; 
right axis). The upward trend means that the low energy com- 
ponent leads the high energies at these long timescales (b) high 
frequency lag-energy spectrum for i/ = [0.96 — 2.98] X 10~^ Hz. 
The nearly zero-lag at low energies, and again at ~ 6.5 keV mean 
that the signal at these energies respond at the same time as the 
reference band. 



The top panel of Fig. [2] shows the energy dependence 
of the lag for the low frequencies where we see the hard lag. 
The lag-energy spectrum is shown in two frequency ranges: 
for the 'low frequency range' {u = [1.65. — 5.44] x lO^** Hz), 
and 'very-low frequency range ' {u = [0.5-1.65] X lO""* Hz). 
The lag spectrum exhibits a general positive trend with in- 
creasing energy. We see a nearly constant zero-lag from bins 
in the 0.3 - 0.9 keV range. At 1 keV, the lag jumps up hun- 
dreds of seconds, meaning that the signal from 1-2 keV is 
occurring hundreds of seconds after the reference band. The 
lag continues to increase with energy up to 10 keV. 

The low frequency and the very-low frequency lag spec- 
tra have the same general shape, but the lag timescales are 
twice as long for the very-low frequency range. For the case 
of propagating viscous fluctuations, this corresponds to vari- 
ability on long timescales occurring further away from the 
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Frequency (Hz) 


Photon Index, F 


u = [0.5 - 1.65] X 10-* 

V = [1.65 - 5.44] X 10"* 

V = [0.96 - 2.98] X 10-^ 

time-integrated 


3.23 ±0.04 
3.12 ± 0.05 
2.89 ± 0.16 
3.228 ±0.009 



Table 2. Power law photon indiees for the eovariance spectra 
in three frequency ranges, and the time-integrated spectrum for 
comparison. 

centre. We also notice that the profile differs at low energies 
between the two frequency ranges. The low frequency lag 
shows a clear step up to ~ 100 seconds, where the very- 
low frequency lag shows a more gradual increase up to the 
constant lag at 1 keV. This suggests that individual com- 
ponents can contribute more at specific frequencies. In this 
case, the soft excess (due to thermal or reflected emission) 
dominates more at the low frequencies, closer to the central 
region, than the very- low frequencies, from further out. 

The bottom panel of Fig. [2] shows the energy depen- 
dence of the lag at higher frequencies from [0.96 — 2.98] x 
10"'' Hz (i.e. the soft lag frequencies). This lag-spectrum has 
a general negative trend, starting at around zero lag in the 
bins from 0.3 - 0.7 keV, and decreasing to ~ —60 seconds 
in the 3-4 keV bin, after which it peaks to zero-lag at the 
energy of the Fe Ka line at ~ 6 keV. Simulations of lags 
between coherent light curves with the same noise level to 
that of the 7-10 keV bin suggest despite low signal-to- 
noise, the lag can be recovered. Though we have taken the 
steps to remove contamination by the Cu-K line of the elec- 
tronics in the 7-10 keV bin, some residuals might effect 
measurements in this bin. 

It is important to note that the high-frequency lag spec- 
trum looks very different than the one at lower frequencies, 
suggesting that difTerent physical processes contribute to the 
hard versus soft lags. 

The high and low frequency lag spectra found 
here with the 1.3 Ms sampl e are similar to those in 
IZoghbi. Uttlev. fc FabianI l|201ll ), but are much more clearly 
defined. Of special interest is the high-frequency feature at 
~ 6 keV, the energy of the iron Ka line. A clear lag sig- 
nature from the Fe L line is not seen in the data, and this 
is likely because of th e thermal reprocessed emission p romi- 
nent at low energies IZoghbi. Uttlev. fc FabiaiilbOlll . and 
see Section [S}. 

3.3 Covariance Spectrum 

The covariance spectrum is a measure of the absolute am- 
plitude of cor related variations in count rate as a functio n 
of energy (see IWilkinson fc Uttlevl[2009l : lUttlev et alll201 j ) . 
Using this approach, we can examine the time-varying por- 
tion of the energy spectrum that is correlated with the ref- 
erence band at a particular frequency range. Note that the 
covariance spectrum is not simply measuring the fractional 
variability amplitude, as in an rms spectrum, but rather is 
proportional to the photon count rate, allowing us to make 
a direct comparison with the time-integrated spectrum. 

Fig. 13] shows the covariance spectrum at the two low 
frequency ranges and the high frequency range. We fit a 
simple power law model at 1 - 3 keV for each spectrum 
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Figure 3. Covariance spectrum for very-low frequencies v = 
[0.5 — 1.65] X 10~^ Hz (blue circle), low frequencies v = [1.65 — 
5.44] X 10"** Hz (red square), and high frequencies u = [0.96 — 
2.98] X 10~^ Hz (black diamonds). At low energies the error bars 
are much smaller than the symbols. We fit a simple power law 
model at 1 - 3 keV. For very-low and low frequencies F ~ 3.2, 
and for high frequencies F ~ 2.9. 

in order to measure the power-law continuum without any 
contamination from the red wing of the Fe Ka line. The 
photon indices are shown in Table [2] We found that the 
photon index hardens with increasing frequency from F = 
3.23 to r = 2.89. We will discuss further implications of this 
result in Section[5l The same spectral hardening is suggested 
when fitting the data from 1.3 - 3 keV alone, though the 
change is not significant. 

Of the three frequency ranges shown here, the high- 
frequency covariance spectrum shows the most prominent 
peak at 0.9 keV, the energy of the Fe L line. The spectrum 
also appears to increase again at 3 - 4 keV, the energy of the 
red wing of the Fe Ka line. These observations further hint 
that the high frequencies are probing the inner accretion 
disc, where relativistic reflection is greatest. 

3.4 Coherence 

Coherence is a measure of the fraction of the variance of one 
process at a given frequency that is a linear transform of the 
other. In other words, a unity coherence shows that one time 
series can completely predict the other. In order to measure 
the lag between two signals, they must have some degree of 
coherence. We c alculate the coherence and error bars using 
E quation 8 from jVaughan fc Nowa^ (|l997f ). Earlier studies 
bv lZoghbi et al.l (|2010l ) examined the coherence as a function 
of frequency between the hard (1-4 keV) and soft (0.3 - 
1 keV) band, and found that for these energies, the coherence 
is high, above ~ 0.75, up to ~ 5 x 10"'^ Hz, where it drops 
sharply due to dominating Poisson noise. 

Fig. [4] shows the coherence in each energy bin rela- 
tive to the reference band (0.3 - 0.8 keV). If the reference 
band is chosen to be the entire energy range (0.3 - 10 keV), 
the shape of the coherence spectrum does not change, but 
the amplitude of the coherence become slightly higher. This 
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Figure 4. The coherence in each energy bin at very-low frequen- 
cies V = [0.5 — 1.65] X 10~* Hz in blue circles, low frequencies 
V = [1.65 — 5.44] X 10~* Hz in red squares, and high frequencies 
u = [0.96 - 2.98] X 10-3 Hz in black diamonds. In Section g] we 
show that the lag can still be accurately recovered despite low 
coherence. 

shows that any structure seen in the coherence spectra is 
not an artefact of the chosen reference band. 

For all of the examined frequency ranges, the coherence 
is 1 below 1 keV, and jumps down at 1 keV, suggest- 
ing that we are seeing contributions from distinct variable 
components. 

In the highest energy bins, the coherence drops, espe- 
cially for high frequencies. In the following sections, we de- 
scribe our Monte Carlo simulations, which show that dilu- 
tions from other uncorrelated components can decrease the 
coherence, but the lag is still recovered. We therefore con- 
clude that while some degree of coherence is necessary, unity 
coherence is not a strict requirement for an accurate mea- 
surement of the lag. 



4 SIMULATIONS 

We conducted simulations to understand the nuances of the 
lag calculation, and to gain clearer insights into the possible 
reasons for the structure seen in the lag spectra. In this 
section, we examine the robustness of the lag measurement 
using Monte Carlo simulations. 

Simulated light curves wore made from rod-noise power 
spectra (PSD) using the Timmcr & Koonig (1995) method, 
which allows for randomness in both the phase and am- 
plitude of the Fourier transform. The underlying red-noise 
PSD has a spectral index of 1.6, consistent with the 
PSD of the data from 0.3 - 10 keV. In order to ac- 
count for possible red-noise leakage into higher frequen- 
cies, the light curves were made ~ 10 times too long, 
and cut do wn to the average light curve length of the ob- 
servations l|Uttlev. McHardv. fc Papadakid [20021 ). The hght 
curves were then scaled to the average count rate and vari- 
ance of all the observations. To simulate a lag, light curve 
pairs were generated from the same random realization of 



the Fourier transform. A frequency-independent time delay 
was then added to the Fourier transform for one light curve. 
The lag and coherence was measured in the high frequency 
range u = [0.96 - 2.98] x 10"^ Hz. 

With these simulated light curves, we test the effects of 
Poisson noise and dilution on the measurements of coherence 
and lag. 



4.1 The effect of Poisson noise 

In order to understand the effect of Poisson noise on the 
measurement of the lag, we simulated 1000 coherent light 
curves pairs that were subjected to increasing amounts of 
noise. A constant time lag of 50 seconds was added to one 
of the light curves in the pair. We calculate the lag for each 
of the 1000 light curve pairs made in the simulation. The 
mean lag was determined from the 1000 measured lags, and 
the standard deviation of the lag distribution calculated. Re- 
gardless of the Poisson noise, the mean lag recovered remains 
consistent with the input, though increasing the noise sig- 
nificantly increases the standard deviation of the recovered 
lags. We test the recovery of the lag errors for 14 simulated 
light curve pairs (i.e. the same number of observations we 
use to calculate the lag), and find that the sizes of these 
error bars are consistent with the data. 

Poisson noise has the effect of adding a component with 
a random phase and amplitude to the Fourier transform of 
our signal. Since the lag is a measurement of the phase lag 
between coherent signals, an incoherent random phase does 
not effect the value of the lag. 

These simulations show that Poisson noise does not ef- 
fect the value of the measured lag as a consistent estimator 
of the intrinsic lag. We find that at high frequencies and 
at high energies where there is low signal-to-noise, it is still 
possible to measure a meaningful lag, as long as appropriate 
error calculations are made. 



4.2 The effect of dilution 

Dilution is another possible effect on the measurement of 
the lag and coherence. If the signals are contaminated by 
an uncorrelated signal or by another correlated light curve, 
the lag measurement will be affected. We check the effect of 
dilution with three scenarios shown in Fig. [5] 

(i) Adding an incoherent component: We simulate 1000 
pairs of coherent light curves with a 50 second time delay. We 
contaminated a fraction of one of the coherent light curves 
with an uncorrelated light curve (different realization of the 
random Fourier amplitude and phase) , and measured the lag 
and the coherence. We find that the coherence drops, but the 
value calculated for the lag is still accurately recovered. 

(ii) Mixing coherent components: We simulated 1000 co- 
herent light curves with a 50 second time delay, and mix one 
light curve with a fraction of the other light curve. We find 
the lag to be a simple average of the mixed components, 
while the coherence remains at 1. 

(iii) Two varying components: We simulate one set of 
1000 coherent light curve pairs offset by 50 seconds, and 
another set of 1000 light curve pairs offset by 20 seconds. 
The light curves offset by 50 seconds do not have the same 
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Figure 5. The recovery of the lag and coherence for diluted signals. When the fraction on the x-axis is 1, the signal is completely 
undiluted, and the 50 second lag is recovered with unity coherence, (i) The light curve with a 50 second delay is contaminated by a 
fraction of uncorrelated signal, (ii) The lag is measured between light curve pairs with 50 second offset that are being mixed together. 
The recovered lag is a simple average of the lags of the different spectral components, (iii) The lag is computed between two variable, 
but uncorrelated light curve pairs, one pair is offset by 50 seconds, and the other is offset by 20 seconds. The coherence decreases with 
more mixing. Shaded regions represent the standard deviation of the 1000 simulated lags. 



realization of the Fourier amplitude and phase as the 20- 
second set. Thus, the two components are variable, but not 
coherent. We dilute the 50-second pairs with an increasing 
amount of the 20-second pairs and find that the recovered 
lag is a weighted contribution from both components. This is 
representative of a scenario in which there are two variable 
components that vary incoherently with each other, have 
different lags and different spectral shapes. This simulation 
shows that the measured lag can be a weighted contribution 
from mixed spectral components, and particularly, that low 
coherence may be an indication of this mixing. 

For clarity. Fig. [5] shows the effect of dilution without 
the addition of Poisson noise into the signal. Even with the 
addition of Poisson noise (equivalent to that seen in the 
data), the effects of dilution are the same. 

These three scenarios show the signal can be recovered 
despite contamination by an incoherent signal. Also, if sig- 
nals are mixed, we will see a proportional drop in the mea- 
surement of the lag difference. We also show that requiring 
a coherence of 1 for a reliable lag measurement is overly 
conservative. 

In the context of the lag spectra for 1H0707-495, we 
see from these simulations that in each energy band, the lag 
is not just due to the reflection or the power law alone, but 
rather from a cont ribution of both. The effect of dilu tion has 
been discussed by IZoghbi. Uttlev. fc Fabiaiil (|201lh , and in 
the following section, we take this idea further by modelling 
the lag spectrum. 



5 SIMULATING THE LAG SPECTRUM 

In the previous section, we showed the effect of dilution on 
the lag spectrum. In this section, we show how the lag spec- 
trum can be decomposed via simulations by accounting for 
dilution in the signal. We offer these simulations as an il- 



lustration of a consistent reffection model that follows the 
data, and not as a statistical fit. We first review our simple 
reffection model for the time-integrated spectrum, and then 
show how the lag-energy spectrum can be modelled in light 
of a refiection scenario. 



5.1 Modelling the time-integrated spectrum 

In this section, we briefly describe our model for the time- 
integrated energy spectrum of the entire 1.3 Ms data set. 
While the energy spectrum is not the focus of this paper, 
we will rely on the spectral model for simulating the lag- 
energy spectrum in the next section. 

Because we are working with such a large amount of 
data, we flnd that our model flt is completely dictated by the 
spectrum below 1 keV, where the signal-to-noise is highest. 
This is ameliorated when including a 2 per cent systematic 
uncertainty to all bins, which is a rough, but low-end, esti- 
mate of the calibration uncertainty. This allowed the high 
energy component (and in particular the clean profile of the 
Fe Ka line) to contribute to the model fit. We also split 
the 1.3 Ms data into 4 ffux-resolved spectra, and fit the 
same models, to understand how the spectra behave in dif- 
ferent states. We find no significant difference in the spectral 
shape aside from changes in the nomalization. While ion- 
ization changes could be expected in different states, so too 
could we expect density change if the accretion rate changes. 
Therefore it is perhaps not surprising that the spectral shape 
remains significantly unchanged between flux states. 

W e fir st follow prev i ous sp ectral work by Fabian et al.l 
l|2012l ) and|D auser et al.l l|2012l ) in fitting the spectrum with 
2 reflection models that are irradiated by a single, steep 
power law. The reflection models cons ist of ionized reflectio n 
from an accretion disc (reflionx ; IRoss fc Fabian|[2005l ). 
which is relativistically smeared (kdblur2f). We also in- 
clude a blackbody to account for both direct disc photons. 
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and thermally reprocessed photons from an irradiated disc. 
Thi s model yields vVdo f = 1940/1528 = 1.27. 

iDauser et al.l 1 20121 ) showed that adding absorption 
from a highly ionized outflow (swindI : iGierlihski fc Pond 
12004 ) significantly improves the fit. In our modelling of 
the entire data set, we also find this to be a good fit 
(x^/dof = 1615/1524 = 1.06), and agree with the possibil- 
ity of a highly ionized outflow. However, as this additional 
component does not significantly change any of the other 
fit parameters, it does not affect our simulations of the lag- 
energy spectrum, and we do not consider it further in this 
paper. 

Given our results from the covariance spectrum dis- 
cussed in Section 13.31 we examined a model in which the 
2 reflectors are irradiated by 2 power laws with different 
spectral indices. The spectral hardening of the covariance 
spectrum at high frequencies suggests that there could be 
two coronal regions, one that originates close to the centre 
and produces a hard power law with short timescale vari- 
ability, and another more extended region which produces 
a steeper power law with longer variations. Both power law 
components likely contribute to the reflection. We simplify 
this complex model of two power laws and two reflectors 
by assuming the same emissivity profile for both reflectors. 
We assume a maximally spinning black hole with a steep 
emissivity that flattens to g = 3 at a break radius of 2.5 rg 
(where rg = GM/c?). This model results in the best fit of 
xVdof = 1638/1526 = 1.07. 

This picture o f an extended coron a is congruous with 
observations from iFabian et al] l|2012l '). when 1H0707-495 
went into a low state, which was interpreted as emission 
produced by an irradiating source within 1 rg of the event 
horizon. The photon index for this low state observation 
was r = 2.7llQjy. In this study, the hard power law in our 
time integrated spectrum is F ~ 2.65, and in the covariance 
spectrum at high frequencies is F = 2.89 ± 0.16. 

Fig. [S] shows our best-fit model with 2 reflectors and 2 
powerlaws. The spectral parameters are shown in Table O 
This model will be used in the following section for simulat- 
ing the lag-energy spectra. 



5.2 Modelling the lag-energy spectrum 

Similar to the Monte Carlo simulations in Section [l] we sim- 
ulate 1000 sets of three identical (but time-shifted) light 
curves to represent three spectral components: a blackbody, 
power law, and reflector. The light curves are simulated 
without Poisson noise, but as noted in Section [4.21 the gen- 
eral effects of dilution on the lag are the same. The signal in 
each energy band is composed of a fractional contribution 
from each of these three components, weighted by their con- 
tributions in the time-integrated energy spectrum shown in 
Fig.H 

We note that the time-integrated spectrum is not the 
ideal measure for the weighting because each of the spec- 
tral components likely contributes more in certain frequency 
ranges than in others. However, we flnd that we can still 
model the lag by making some assumptions about the vari- 
ability of the components in the time-integrated spectrum: 

(i) The blackbody component is likely composed of both 
direct disc photons (that would be appreciable in the low 




Energy (keV) 

Figure 6. Time-integrated energy spectrum for the entire 1.3 
Ms exposure. The data are modelled with a blackbody (green 
dot-dash), 2 power laws (thick blue and thin red sohd) and 2 
reflection spectra (low ionization reflector, J ~ 50 irradiated by 
the soft power law as a thin red dotted line; and high ionization, 
5 ^ 500, irradiated by the hard power law as the thick blue dotted 
line). We use the fraction of each component to the total model 
in order to weight the light curves in our simulations of the lag- 
energy spectrum. 



Component 


Parameter 




Value 


Galactic absorption 
intrinsic absorption 


AfH(1022 cm- 
AfH(1022 cm- 


2^ a 


0.058 

Qo-, +0.005 
u.uoi „Q 002 


blackbody 


Abb X 10-5 
fcTbb (eV) 




4 5 +^-^ 
-0.6 

88.812-2 


power law 1 


X 10-* 

Fi 




2.3 

2 66+° °* 


power law 2 


Ay^ X 10-* 

r2 




7 4 +0-4 
-1.3 

3 30+° ° 
■^•■^^-0.02 


kdblur2f 


'nn(rg) 
rout(rg) " 

in 

^break ('"g ) 
92 ° 




1.235l°:°i 
10.0 1°, 
400.0 
46.8 ±0.4 
2 49+0-03 

^•^="-0.02 

3.0 


reflionx 1 


A.efli X 10-6 

4efli (erg cm 


s-1) 


0.13 ±0.02 

484 1°, 
13.4 l°i 


reflionx 2 


A,efl2 X 10-'' 

Crefla (erg cm 


s-1) 


43.0 +11 

AQ 7 +0-4 


xVdof 


1638/1526 = 


1.07 





" frozen 

Table 3. Best-fit spectral parameters for the time-integrated 
spectrum. 
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Low Freq. Lags (s) 


High Freq. Lags (s) 


blackbody 








Thard 


1000 


-90 


Tsoft 


400 




reflection 




-10 



Table 4. Selected time delays for each component in the high 
and low frequency lag simulations. The smaller the value of the 
lag, the earlier it responds. 



frequency signal) and thermal photons that are produced 
by the irradiating continuum heating up the surface of the 
disc (discernible in the high frequency signal). Because we 
cannot distinguish between the low and high frequency con- 
tributions in the time-integrated spectrum, we use the entire 
blackbody when calculating the weighting for both the high 
and low frequency lags. 

(ii) We have seen in Section [3.31 that the power law does 
not contribute the same at all frequencies, and there is likely 
mixing between at least two power-law components at dif- 
ferent frequencies. For the low frequencies, the fluctuations 
propagate inwards and are upscattered in the soft power 
law emitting corona first, and then later in the inner hard 
corona. In our simulations of the low frequency lag spec- 
trum, we include the signals from both the hard and soft 
power-law components. For the high frequencies, where we 
are presumably probing closer to the inner regions, the hard 
power-law component dominates, and therefore in our high 
frequency model, we only simulate the hard power-law com- 
ponent. 

(iii) We assume that the lag due to reflection is most ap- 
preciable on short timescales, and so we only include the re- 
flection component when modelling the high frequency lag 
spectrum. It is expected that reflection contributes some- 
what at low frequencies as well, but the signal is dominated 
mostly by the variability from propagating viscous fluctua- 
tions. 

The time delay between the contributing spectral com- 
ponents is assumed to be constant for this simple model 
(i.e. assuming a delta function transfer function). We input 
reasonable expectation values for the time delays which are 
obtained from the geometry implied by our model. Tabled] 
gives the time delays for both the low and high frequency lag 
spectral models. We reiterate that the smaller the value of 
the lag, the earlier it responded. For instance, at low frequen- 
cies, the power law lags behind the blackbody, and therefore 
its input lag is greater. At high frequencies, the power law 
leads the reflection, and so its input lag is more negative than 
the reflection lag. In other words, at high frequencies, the 
power law responds flrst, followed by reflection 80 seconds 
later, and finally followed by the reprocessed blackbody 10 
seconds after the reflection. Fig.Qillustrates how the signals 
from each component are weighted by the time-integrated 
spectrum to produce the simulated lags. 

The time lag is then calculated in the usual way, with 
respect to the reference signal from 0.3 - 0.8 keV. The top 
panel in Fig. [8] shows the simulated lag spectrum in the 
frequency range v = [0.5 - 1.65] X IQ-" Hz. Notice how the 
step up from seconds to ~ 200 seconds arises naturally due 
to the time delays between the thermal emission in the disc 
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Figure 7. Illustration demonstrating how the high frequency lag 
spectrum is simulated. (i) The input time delays for the hard 
power law (green), reflection (red) and blackbody (blue) as a 
function of energy. The power law responds first, followed by the 
reflection 80 seconds later, and the reprocessed black body 10 sec- 
onds after that, (ii) The fraction of the hard power law, reflection 
and blackbody to the total time-integrated spectrum, (iii) The 
simulated lag that is compared to the data in bottom panel of 
Fig. [8] Each signal in (i) is weighted by its contribution in (ii), 
and the total time lag is then summed to make the simulated lag 
spectrum in (iii). 



and the continuum emission in the corona. This can be easily 
explained by the fractional contribution of the blackbody 
and two power-law components in each band compared to 
the reference. For instance, at 3 keV, the signal is roughly 
composed of two parts soft power law (with lag = 400 s) 
and one part hard power law (lag = 1000 s), amounting 
to an average lag of ~ 600 seconds. This lag is compared 
to the reference band, which is roughly composed of 30% 
blackbody (lag = s), 10% hard power law and 50% soft 
power law, amounting to an average lag of ~ 300 s. The lag 
at 3 keV, is therefore calculated to be ~ 300 s. 

The bottom panel of Fig. [S] shows the lag in the fre- 
quency range = [0.96-2.98] X 10"^ Hz. The model follows 
the general shape of the data, and suggests that reflection 
is playing a major role in diluting the intrinsic lag of the 
power law. However, the model is particularly unsatisfac- 
tory in the 3-4 keV bin. At 3-4 keV, the red wing of the 
iron Kq line dominates the reflection spectrum. This is the 
reflected emission from the innermost regions of the disc, 
and therefore, it is the signal that responds first. To account 
for this, we input a time lag of -75 s for the refiection in 
this bin, which amounts to the reflection lagging behind the 
power law by just 15 seconds. Therefore, the effect of dilu- 
tion on the power law due to reflection is less in this band, 
making the lag appear to be lower. Fig. [9] shows the results 
of the simulations with the addition of a fast response from 
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Figure 8. Monte Carlo simulations of the low and high frequency 
lag-energy spectrum (black dashed lines) where the weighting of 
the signal from the blackbody, power law and reflection compo- 
nents are determined by their fractional contribution to the time- 
integrated energy spectrum. Time delays between all components 
are constant. 



the red wing of the Fe Ka line. The addition of the response 
from the red wing significantly improves the fit to the data. 



6 DISCUSSION 

6.1 The Low Frequency Lags 

The large positive lag at low frequencies has been explained 
as the vi scous propagation of accretion rate fluctuations in 
the disc (lArevalo fc UttlevlbOOd : I Wilkinson fc UttlevllioogI : 
lUttlev et al.ll201ll ). If the energy released from accretion is 
dissipated in the corona (e.g. through magnetic fleld lines), 
fluctuations in the disc could cause correlated variability in 
the corona. Our analysis of the hard lag spectrum is con- 
sistent with this model. We find zero lags with respect to 
the reference band (0.3 - 0.8 keV) at energies below 1 keV 
(where the blackbody is appreciable), and above 1 keV, 
where the power law begins to dominate, we compute a lag 
behind the reference of 200 s up to 600 s. Considering the ef- 
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Figure 9. Monte Carlo simulations for the high frequency lag 
spectrum. This figure has the same model as the bottom panel 
of Fig. [8] however, we now account for a faster response from the 
red wing of the Fe Ka line. In our model, the red wing of the line 
responds 15 seconds after the continuum emission, while the bulk 
of the reflected emission responds 80 seconds after the continuum. 



feet of dilution in our simulations, we find that the intrinsic 
time delays are likely to be greater (~ 400 s to 1000 s). 

Furthermore, our frequency-resolved analysis of the low 
frequency lag-energy spectrum (top panel. Fig. [2]), shows 
that for variations on longer timescales ('very-low' frequen- 
cies, v — [0.5 - 1.65] X 10"* Hz), the blackbody compo- 
nent contributes less. These very-low frequency variations 
presumably originate from further regions from the centre, 
where the temperature of the disc is cooler, and may be 
out of the spectral sensitivity of XMM-Newton. At higher 
frequencies, {y = [1.65 — 5.44] x 10~* Hz), the variations 
originate closer in, where the disc is hotter, and therefore, 
the contribution from the blackbody is greater. Our simula- 
tions also confirm that if we increase the fractional contri- 
bution of the blackbody (as we expect for smaller radii in 
a Shakura-Sunyaev disc), we describe the general shape of 
the lag-energy spectrum for the 'low' frequencies in the top 
panel Fig. H 

The low frequency lag spectrum shows a steady increase 
above 1 keV. This upward trend has b een observed in ste l- 
lar mass black holes, such as Cyg X -1 jNowak et al.lll999l ). 
iKotov. Churazov. fc Gilfanovl (j200ll ) argued that these lags 
could be due to radial dependence of the spectral index, in 
which it is hardest in the innermost regions, and becomes 
softer at further distances. Perturbations in the accretion 
flow that are introduced at different radii (and therefore dif- 
ferent timescales) propagate inward on diffusion timescales, 
which cause the observed power law shape of the lag spec- 
trum. This model for Galactic black holes can also explain 
the upward trend seen in the lag spectrum of 1H0707-495. 
The data are well described by our simulations, which fol- 
low the Kotov model of the soft power law responding first 
followed by the hard power law from the innermost regions. 

Finally, this interpretation of an extended corona is con- 
sistent with relativistic models by Wilkin s fc Fabian ( 201 J ) 
that show that a majority of the power law comes from a 
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small region within 5 rg, while the rest comes from further 
out. 

6.2 The High Frequency Lags 

The negative lags at high frequencies have been explained as 
the light travel time between continuum emission produced 
in the corona and reflected emission from the inner regions 
of the disc. The soft lag spectrum in bottom panel of Fig.[5]3 
is generally consistent with this model in that variations at 
energies dominated by the reflection spectrum are tempo- 
rally coincident with those of the reference band, while the 
energies dominated by the power-law component lead those 
variations. This result is consistent with previous results, 
and now with the new long dataset, we are able to define 
the lag spectrum much more clearly, revealing structure at 
the Fe Ka line. 

An important result from the simulations in Fig. [S] is 
to show that the actual time delays between the power law 
emission and the reflection are not necessarily the same as 
the lags we measure because of the effect of dilution. We find 
that for much of the refiection, the actual light travel lag is 
more on the order of 80 seconds, corresponding to a coronal 
height of '-^ 4 r-g (assuming a black hole mass 2 x 10^ Mq). 
At 3 - 4 keV, the red wing of the iron Ka line, we find 
that the light travel time is ~ 15 seconds, corresponding to 
emission from within 1 rg of the event horizon. Future work 
using realistic transfer functions, which take into account 
the spatial extent of the source and the inclination of the 
disc, will improve our soft lag spectral models, and give us 
a better understanding of the geometry of the innermost 
regions closest to the black hole. 



7 CONCLUSIONS 

We have presented the timing analysis of 1H0707-495 using 
1.3 Ms of data taken over 10 years with XMM-Newton. We 
have simulated the lag spectra using a model that accounts 
for dilution from different spectral components in calculat- 
ing the lag. This simple model for the hard and soft lag 
spectra shows complicated features that can be explained 
by reflection. We have also shown through an analysis of 
the covariance spectrum and the time-integrated spectrum, 
that there is evidence for a more complicated continuum 
emitting region, which produces a harder power law closer 
to the centre, and a softer power law at larger radii. Further 
work needs to be done to understand the complexities of the 
power law emitting region. This, and future work with more 
realistic transfer functions for the response of the accretion 
disc will be beneflcial in uncovering the precise origin of the 
X-ray emission in the local environments of supermassive 
black holes. 
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APPENDIX : CHECKING THE BACKGROUND 
FOR LAGS 

To ensure that the lags in Fig. [1] are intrinsic to the source, 
we replicated the analysis on background light curves, and 
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Figure A. Background lag-frequency spectrum in black, with 
source lag-frequency from Fig. [T] in grey for reference. While the 
general shape of the source and background lag-frequency spectra 
are similar, the background lags are not significant. We perform 
this background analysis to confirm that the lags in Fig. [l] are 
intrinsic to the source. 

found no significant lag. Tiie same Fourier analysis, as de- 
scribed in the previous section, was performed using the 
hard and soft background light curves of the 14 observations. 
The results are shown in Fig. [X] The general shape of the 
background lag-frequency spectrum is similar to that of the 
source, but the amplitude of the lags is significantly smaller. 
This could perhaps be caused by contamination of the back- 
ground by source photons. A complete study of these mild 
background fluctuations is beyond the scope of this paper, 
but we do conclude that the lags seen in Fig. [1] are intrinsic 
to the source. 
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